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Abstract 

^ • The pericenter advance of a test body by a central mass is one of the classical tests of General Relativity. 

^ ■ To day this effect is measured with radar ranging by the perihelion shift of Mercury and other planets, in 

t^X}, the gravitational field of Sun, with a relative accuracy of the order of 10^^-10^''. In this paper we explore 

the possibility of a measurement of the pericenter advance in the gravitational field of Earth by analyzing the 
laser-ranged data of some orbiting, or proposed, laser-ranged geodetic satellites. Such a measurement of the 



X 



^ ' perigee advance would place limits on hypothetical, very weak, Yukawa-type components of the gravitational 



interaction with a finite range of the order of 10^ km. Thus, we show that, at the present level of knowledge of the 
orbital perturbations, the relative accuracy, achievable with suitably combined orbital elements of LAGEOS and 
LAGEOS II, is of the order of 10"'^. With the corresponding measured value of ^^^3 ~^ , by using 77 = 4/3 — 7 — 3 
from Lunar Laser Ranging, we could get an estimate of the PPN parameters 7 and /3 with an accuracy of the 
order of 10^^ — 10^^. Nevertheless, these accuracies would be substantially improved in the near future with 
the new Earth gravity field models by the GHAMP and GRAGE missions. The use of the perigee of LARES 
(LAser RElativity Satellite), with a suitable combination of orbital residuals including also the node and the 
perigee of LAGEOS II, would also further improve the accuracy of the proposed measurement. 
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1 Introduction 

The well known secular shift of the pericenter a; of a test body induced by the Schwarzschild 
metric generated by a static, spherically symmetric distribution of mass [Misner et al, 1973; 
Ciufolini and Wheeler, 1995] is one of the classical tests of Einstein's General Relativity. 

1.1 The Solar System tests 

The advance of pericenter, which can be expressed in terms of ^+^^~^ ^ where f3 and 7 are 
the Eddington-Robertson-Schiff PPN parameters [Will, 1993], has been measured in the Solar 
System, on the Mercury perihelion, by radar signals transmitted from Earth to Mercury and 
back to Earth [Shapiro et al, 1972], yielding 2±^^ = 1.005 ± 7 x 10"^ The contribution of 
possible systematic errors raises the uncertainty to about 2 x 10^^. From the analysis of ten 
years of data the relative error published in [Shapiro et al, 1976] amounts to 5xl0~^. The 
major sources of systematic error in this measurement are the poorly-known variations in the 
topography of the planet's surface and the uncertainties in the radar scattering law [Shapiro, 
1990; Pitjeva, 1993]. For a review of periastron advance measurements on binary pulsars see 
[Will, 2001] and references therein. 

Several space missions have been proposed in order to measure, among other effects, this 
phenomenon by using a variety of techniques: the most recent and promising are SORT, IPLR 
and ASTROD [Ni, 2001]. 

More accurate measurements of might be performed in the future by means of the 

ESA BepiColombo Mercury orbiter [Balogh et al, 2000; Milani et al, 2001a; Milani et al, 
2001b]. In [Milani et al, 2001b] it is claimed that the BepiColombo mission, which should be 
launched in 2009, should allow to measure 7 and (3 with a relative accuracy of the order of 

2 X 10"^ 

1.2 The proposed experiment 

In this paper we propose to perform a measurement of the general relativistic "gravitoelectric" 

pericenter advance in the gravitational field of the Earth by using some suitable combinations of 
the orbital residuals of existing or proposed, spherical, passive, geodetic laser-ranged satellites 
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[lorio, 2002a] with particular emphasis to LAGEOS and LAGEOS II in order to exploit the rel- 
evant experience obtained with the gravitomagnetic Lense-Thirring measurements in [Ciufolini 
et al, 1996; 1997; 1998]. In [Ciufolim and Matzner, 1992] a determination of the LAGEOS 
general relativistic gravitoelectric perigee shift in the field of the Earth is reported, but the 
accuracy amounts only to 2 x 10~^. 

This paper is organized as follows: in section 2 we compare the present experimental accu- 
racy in satellite laser-ranging measurements with the relativistic expression of the perigee shift 
of the two LAGEOS satellites. In section 3 we analyze some of the most important sources of 
systematic errors. The accuracy achievable in measuring f3 and 7 is also discussed. By detect- 
ing the relativistic perigee rate of the proposed LARES [Ciufolini, 1998] laser-ranged satellite, 
in conjunction with the node and the perigee of LAGEOS II, it might be possible to further 
improve the accuracy in the measurement of these parameters over a long enough time span. 
This is the subject of section 4 together with an analysis of the potential role of the existing 
laser-ranged satellite Starlette. Section 5 is devoted to the conclusions. 

2 The relativistic perigee precession of LAGEOS type 
satellites 

As known, in the slow-motion and weak-field approximation, the Schwarzschild metric gen- 
erated by a static, spherically symmetric distribution of mass-energy induces an additional 
post-Newtonian "gravitoelectric force" which acts on the orbit of a test body by shifting its 
pericenter; in the PPN formalism the pericenter rate can be written as: 

. _ SnGM 2 + 27-/? 

"^GE — -5—7:; X 

in which G is the Newtonian gravitational constant, c is the speed of light in vacuum, M is the 
mass of the central object, a and e are semimajor axis and eccentricity, respectively, of the orbit 
of the test body and n = \JGM /a^' is its mean motion. In the following we define v = ^^^3"^ ■ 
General Relativity predicts that the perigee shifts for LAGEOS and LAGEOS II amount to 
3,312.35 miUiarcseconds per year (mas/y in the following) and 3,387.46 mas/y respectively. 

The actual experimental precision allows for detecting such rates for both the LAGEOS 
satellites [Ciufolini, 1996]. Indeed, for the perigee the observable quantity is r = eojjj] at 
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present the RMS of the residuals of the best orbital fits amounts to Srg^p ~ 1 cm, or less, for 
the two LAGEOS satellites, over several orbits and with a set of force models (i.e. not including 
modeling errors). Since the LAGEOS eccentricity is ei = 4.5 x 10"'^ the accuracy in detecting 
the perigee is Scul^p ~ 5rexp/ (ejai) ~ 37 mas. So that, over 1 year the relative accuracy in the 
measurement of the relativistic perigee shift would be of the order of 1 x 10~^. For LAGEOS II 
this measurement accuracy is better than for LAGEOS. Indeed, the LAGEOS II eccentricity is 
en = 1.4 X 10~^, thus 5a;"p ~ 12 mas; this may yield an accuracy of the order of 3 x 10~^ over 
1 year. Over 8 years, by assuming the same fit error as before, it would amount to ~ 4 x 10^''. 
These considerations rule out the possibility of directly using the perigee of the other existing, 
spherical, geodetic laser-ranged satellites Etalon-1, Etalon-2, Ajisai, Stella, Westpac-1 because 
their eccentricities are even smaller than that of LAGEOS. On the contrary, Starlette has 
an eccentricity of the order of 2 x 10~^; however, since it orbits at a lower altitude is more 
sensitive than the LAGEOS satellites to atmospheric drag and to Earth's zonal harmonics, so 
that would be difficult to process its data at an acceptable level of accuracy (See also section 
4). Accordingly, in order to detect the gravitoelectric relativistic shift in the gravitational field 
of Earth, we will focus on the perigee of LAGEOS and especially of LAGEOS II . 

3 The systematic errors 

The perigee of an Earth sateUite is a " non-clean" orbital element in the sense that it is affected 
by a large number of aliasing classical forces inducing systematic errors which must be carefully 
identified and analyzed. 

3.1 The static geopotential error 

The most important source of systematic errors in such a measurement is represented by the 
mismodeling of the classical perigee precession induced by the errors in even zonal harmonics 
of the Earth's gravitational field [Kaula, 1966]. This error is really critical because the resulting 
aliasing trend cannot be removed from the data and nothing can be done about it apart from 
assessing the related error as reliably as possible. 

By using the covariance matrix of the EGM96 Earth gravity model [Lemoine et al., 1998] 
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and adding the correlated terms in a root-sum-square fashion up to degree / = 20, we obtain 
for LAGEOS a systematic error 5z//i/zonais = 8-1 10^^ whereas for LAGEOS II we have 
5^/^zonais = 1-5 X 10~^. Siuce the major source of error lies in the uncertainty of the first two 
mismodeled even zonal harmonics 5J2 and 5J4, in order to eliminate most of the static and 
dynamical even zonal terms of the geopotential, following the method in [Ciufolini, 1996] for the 
measurement of the gravitomagnetic Lense-Thirring effect, we search for suitable combinations 
of the orbital residuals of the existing laser-ranged satellites. In Tab.|l] we report the most 
promising combinations: their general form is 

N 

5lu^^ + ^ Ci^fi' -I- c^^_^_lSuJ^ = xgr z^, (2) 

i=l 

in where is the number of the nodes of different laser-ranged satellites employed, xqh is the 
slope, in mas/y, of the relativistic trend to be measured, in Tab.|l] Su/uzonais is the systematic 
error induced by the even zonal harmonics up to degree / = 20 calculated with the EGM96 
covariance matrix. The coefficients Cj are determined in terms of the orbital parameters of the 
satellites entering the combinations. 



Table 1: Combined residuals: numerical values 











J 








Cl 




C3 


C4 


xgr (mas/y) 


S^/ Z^zonals 


1 


-0.87 


-2.86 








3,387.46 


6.59 X 10-3 


2 


-2.51 


-4.37 





2.51 


11,704.92 


1.1 X 10-3 


3 


-1.96 


-3.69 


0.037 


1.37 


7,928.51 


8.1 X 10-^ 



It is important to stress that the use of the LAGEOS satellites, due to their altitude, makes 
our measurement substantially insensitive to the errors in the zonal harmonics of degree / > 20, 
so that our estimates of Su/u^onais presented here are valid even in the case that the EGM96 
covariance matrix for degrees higher than / > 20 would not be accurate enough. 

In addition to LAGEOS and LAGEOS II we have only considered Ajisai [lorio, 2002a] 
since it is well tracked, contrary for example to the Etalon satellites, and it would be less 
demanding than for the other satellites to reduce its laser-ranged data to a level of accuracy 
comparable to that of LAGEOS and LAGEOS II. Moreover, the other laser-ranged satellites 
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orbit at lower altitudes, therefore they are more sensitive to terms of the geopotential of higher 
degree. Consequently, as confirmed by numerical calculations, the inclusion of their data in the 
combined residuals would increase 5z//z/zonaisB 

It is important to point out that the uncertainties quoted in Tab.|l] for 5v/vzon&\s will be 
reduced in the near future when the new gravity models from the CHAMP and GRACE missions 
will be released [Pavlis, 2002a]. 

3.2 The time— dependent systematic errors 

In regard to the evaluation of the impact of the other sources of systematic errors, the role of 
the coefficients q entering the combinations has to be taken into account, indeed, according to 
their magnitude, they can reduce or emphasize the effects of the perturbing forces acting on 
the orbital elements weighted by them. 

For example, for the combination 2 of Table 1, over 1 year, the impact of the error in 
measuring the perigee rate of LAGEOS amounts to about 2.5 x (3.1 x 10"'^) = 7.7 x 10^'^ while 
for combination 3 it is about 1.37 x (4.7 x 10"^) = 6.4 x 10~^. 

In this paper we only consider in detail the combination 1 of Table 1 

5u^^ - 0.87 X - 2.86 x Sti^ = 3, 387.46 x u, (3) 

so that we exploit the insight acquired with the Lense-Thirring LAGEOS experiment [Ciufolini 
et al, 1996; 1997; 1998]. 

The long-period harmonic perturbations, according to their periods P and to the adopted 
observational time span Tobs, may turn out to be less insidious than the mismodeled secular 
perturbations due to the zonal harmonics of the geopotential since, if P < Tobs and Tobs = 
nP, n = 1, 2, ... they average out; if their periods are shorter than the time span they can be 
viewed as empirically fitted quantities which can be subsequently removed from the signal. 

The results and estimates recently obtained for the LAGEOS gravitomagnetic experiment 
can be applied to our proposed measurement. Indeed, the solid Earth and ocean tidal pertur- 
bations on LAGEOS and LAGEOS II and their impact on the Lense-Thirring effect have been 

^ About combination 3, it should be noted that, in order to obtain more rehable and accurate estimates of 
the systematic error due to the even zonal harmonics of the geopotential, the calculation of the error should 
include terms with degree higher than / = 20 due to the sensitivity of Ajisai to such higher degree terms. 
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studied in [lorio, 2001; lorio and Pavlis, 2001; Pavlis and lorio, 2002], while the role played 
by the non-gravitational perturbations has been analyzed in [Lucchesi, 2001; 2002]. Regarding 
the atmospheric drag, whose impact might be thought of as a serious drawback, especially for 
the perigee of LAGEOS II, in [Ciufolini et ai, 1997] it is shown that it essentially averages out 
over an orbital revolution. 

In Tab.0 we summarize the results obtained with the proposed combination using LAGEOS 
and LAGEOS II for a 8-year time span with orbital fits, arcs, each of 7 days. In assessing 
the total systematic error we have accounted for the fact that the gravitational errors are not 
independent, indeed we have simply summed them up; then we have added these uncertainties 
and the other independent systematic errors in a root-sum-square fashion. In Tab.|^ we can 
observe that the error budget is mainly dominated by the even zonal harmonics and by the 
LAGEOS II perigee measurement error. In regard to the non-gravitational errors, according 
to the results of [Lucchesi, 2001] for the direct solar radiation pressure and the Earth's albedo 
and to the results of [Lucchesi, 2002] for the thermal thrust perturbations and the asymmetric 
reflectivity, the corresponding uncertainty would amount to almost 1 x 10~^ over 7 years. In 
obtaining this result a very pessimistic approach has been adopted by assuming a mismodeling 
of 20% for all the perturbing effects except for the direct solar radiation pressure which has 
been assumed to have an uncertainty at the 0.5% level. However, we stress that only the 
Earth's thermal thrust, or Yarkovski-Rubincam effect, induces a mismodeled linear trend whose 
impact would amount to about 1 x 10~^: the other non-gravitational forces are time-varying 
with known periodicities and can be fitted and removed from the signal as done for the tidal 
perturbations. 

Table 2: Preliminary error budget: Tobs = 8 years. At = 7 days. 

Even zonal harmonics (par. 3.2) 6.59 x 10^'^ 

J3 geopotential [Pavlis and lorio, 2002] 3.2 x lO^'^ 

Tides [lorio, 2001; PavUs and lorio, 2002] 4.4 x 10"^ 

Non-gravitational effects [Lucchesi, 2001; 2002] lxlO~^ 

Measurement error in LAGEOS II perigee (sect. 2) 4 x 10^^ 



Total systematic error 7.3 x 10 



-3 
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Forthcoming improvements in the Earth gravity field modehng achievable with the CHAMP 
and GRACE data have been recently extensively studied and reported in [Pavlis, 2002a]. They 
will yield important improvements in the accuracy of our experiment. We finally observe that 
we already have the required laser-ranging data since, at present, we have almost 10 years of 
LAGEOS and LAGEOS II data. 

3.3 The parameters /3 and 7 

The results obtained for the combination 1 examined here, together with 61] = 8x 10^^ reported 
in [Anderson and Williams, 2001] for the combination 77 = 4/5 — 7 — 3 [Nordtvedt, 1968; 1991], 
would allow to measure /3 = + |i/ + | at the level ~ 3 x 10~^ independently of other 
measurements of 7; this result, which is of the same order of magnitude of that obtained with 
radar ranging [Shapiro, 1990], should be compared with the most recent 513 — 4.7 x 10~^, 
obtained with LLR data [Anderson and Williams, 2001]. 

The parameter 7, that may be written as 7 = i?7 + yZ/ — |, would be measured less precisely: 
^7 = 1.2 X 10^^. However, 7 can be directly measured via light deflection or radar time delay. 
For example, in [Froeschle et al., 1997] the result 7 = 0.997±3 x 10~^ is based on the astrometric 
observations of electromagnetic waves deflection in the visible frequency. By using the Viking 
time delay [Reasenberg et al, 1979], 7 = 1.000 ±2 x 10~^ was obtained; the quoted uncertainty 
allows for possible systematic errors. A more recent measurement based on the time delay with 
the NEAR spacecraft [Elliott et al, 1998] claims an uncertainty < 10~^. 

For an updated overview of the present status of the measurements of the PPN parameters, 
see [Will, 2001] 

4 The role of LARES and Starlette 

The proposed LARES satellite [Ciufolini, 1998] will enable us to perform several general rel- 
ativistic tests and geophysical measurements. In this section we investigate its role for the 
measurement of the gravitoelectric perigee advance that would be well detectable on LARES 
thanks to its proposed relatively large eccentricity: clr — 0.04. 
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A possible observable is the following combination of orbital residuals 

5u^^ + ci5u^^ + ca^fi" = xgkv, (4) 

with 

ci = -4.71, (5) 

C2 = 2.26, (6) 

xgk = -12, 117.4 mas/y. (7) 

The coefficients Cj depend on the orbital parameters of the selected satellites and have been 
calculated corresponding to the proposed orbital parameters of LARES: olr = 12, 270 km, 
zlr = 70 deg and clr = 0.04. Let us now calculate the systematic errors induced by the even 
zonal harmonics of the geopotential, according to the EGM96 gravity model. Eq. would 
allow to cancel out the errors due to the first four even zonal harmonics; the impact of the 
remaining ones, with degree higher than four, amounts to an uncertainty of about 6 x 10~^. 
This error will further reduce when the new gravity models from the CHAMP and GRACE 
missions will be available. An interesting feature of eq. @) is that the impact of the non- 
gravitational perturbations would produce an uncertainty of about 3.9 x 10^'^ over a time 
span of 7 years, according to the very conservative estimates discussed above; in particular the 
Earth thermal thrust, i.e. the so-called Yarkovski-Rubincam effect, would induce a mismodeled 
secular trend with an uncertainty of only 6.4 x 10~^. The inclusion of LARES would thus 
represent an improvement with respect to the LAGEOS and LAGEOS II scenario previously 
outlined, especially when, in the near future, the impact of the mismodeled non-gravitational 
perturbations will increase relatively to the gravitational perturbations thanks to the more 
accurate Earth's gravity field models. 

Among the existing laser-ranged satellites, Starlette, with its relatively large eccentricity of 
0.0204, seems to be a natural candidate for the measurement of the perigee advance. Unfor- 
tunately, it orbits at an altitude much lower than the LAGEOS satellites (astr = 7,331 km) 
and thus the mismodeling of the even zonal harmonics of the geopotential would induce large 
systematic errors in the measurement of the perigee rate. It turns out that the inclusion of the 
perigee of Starlette in the combined residuals if, on the one hand, provides a further observable 
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and thus would allow to cancel out the error of an additional mismodeled even zonal harmonic, 
on the other hand it raises the total error due to the geopotential, even with the LARES data. 
Indeed, in a combination including the nodes of LAGEOS and LAGEOS II and the perigees of 
LAGEOS II and Starlette, the relative error due to the geopotential amounts to only 6.2 x 10~^. 
This rules out the possibihty of using the Starlette data, at present level of knowledge of the 
Earth's gravitational field, to improve the measurement of the relativistic perigee advance. 

5 Conclusions 

The main features of our proposed measurement of perigee advance using laser-ranged satellites 
are the following: 

• It will provide a test of a general relativistic effect, i.e. the pericenter advance, in the 
gravitational field of Earth, complementary to similar tests in the field of Sun. Indeed, 
such a measurement in the gravitational field of Earth would place stronger limits on 

possible, very weak, Yukawa-type components of the gravitational interaction with a 
finite range of the order of 10^ km [lorio, 2002b] . 

• For a suitable combination of the orbital residuals of the nodes of LAGEOS and LAGEOS 

II and perigee of LAGEOS II, the relative error due to the mismodeled classical precessions 
by the static even zonal harmonics of the geopotential, which is the most insidious source 
of systematic error, amounts, according to the EGM96 Earth gravity model, to about 
6 X 10~^. The total error in our proposed combination of orbital elements, also due to 
the other time-dependent systematic errors, amounts, over an observational time span of 
eight years, to almost 7 x 10~^. 

• In the near future the knowledge of the Earth gravitational field will increase thanks to the 
data from the CHAMP and GRACE missions. Moreover, the rms accuracy in measuring 
the position of the LAGEOS satellites should reach the few mm level. Consequently, 
there will be an increase in the accuracy of our measurement. 

• In regard to the terrestrial space environment, which induces a large number of orbital 
perturbations on the satellites, it is very well studied and the action of these perturbations. 
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especially on the LAGEOS satellites, is very well modeled. The data processing would be 
reliably based on the large experience collected in several years of analysis of the LAGEOS 
orbits [Pavlis, 2002b]. Moreover, it would be relatively fast and easy to include the future 
improvements in the force modeling in our analysis, e.g. the improvements due to a better 
experimental and theoretical knowledge of the rotation rate and spin axis orientation of 
LAGEOS type satellites. 

• The proposed measurement, as an additional outcome, from in conjunction with 
77 = 4/3 — 7 — 3 from Lunar Laser Ranging, would allow to obtain the PPN parameters 
7 and P in the gravitational field of Earth (see above) with an accuracy at the level of 
10"^ — 10~^. This accuracy would improve with the new Earth's gravitational models 
from CHAMP and GRACE. 

• Starlette data, although this sateUite has a relatively eccentric orbit, should not be in- 
cluded in the combined residuals. Indeed, Starlette would raise the systematic error in 
the gcopotential due to its lower altitude with respect to the LAGEOS satellites and thus 
due to its higher sensitivity to the even zonal harmonics of the terrestrial gravitational 
field. 

• The use of the perigee of the proposed LARES satellite, for a suitable combination of 
observables including also the node and the perigee of LAGEOS II, would allow to reduce 
the systematic relative error due to the mismodeled even zonal harmonics of the gcopo- 
tential to about 6xl0~^, according to EGM96. This proposed measurement would also 
take advantage of the new Earth gravity models from CHAMP and GRACE. According to 
conservative estimates, the impact of the non-gravitational perturbations would amount 
to about 4 X 10~^ over 7 years; in particular the impact of the mismodeled Earth thermal 
thrust, or Yarkovski-Rubincam effect, would amount to only about 6 x 10~^. 
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